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Q ; Abstract 



Realistic nucleon-nucleon interaction induce correlations to the nuclear many-body system, which 
lead to a fragmentation of the single-particle strength over a wide range of energies and momenta. 



> 

qq \ We address the question of how this fragmentation affects the thermodynamical properties of 

O ' 

I/"} ■ nuclear matter. In particular, we show that the entropy can be computed with the help of a 

O ' 

. spectral function, which can be evaluated in terms of the self-energy obtained in the Self-Consistent 

o 

Green's Function approach. Results for the density and temperature dependences of the entropy 
per particle for symmetric nuclear matter are presented and compared to the results of lowest 

O : 

order finite temperature Brueckner-Hartree-Fock calculations. The effects of correlations on the 

Ch ■ 

calculated entropy are small, if the appropriate quasi-particle approximation is used. The results 

•rH , 

demonstrate the thermodynamical consistency of the self-consistent T-matrix approximation for 

5-H ' 



the evaluation of the Green's functions. 
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I. INTRODUCTION 



The entropy is a key quantity in the study of the thermodynamical (TD) properties of 
fermionic systems. In the context of correlated Fermi liquids, the study of the entropy has 



1 hvii triggered by the experimental and theoretical studies of 3 He [llO, [3 ,13, [3, f|- The spot-ilk- 
leat of this system is experimentally known to have a non-trivial temperature dependence 
?! of the type T 3 lnT. Such a non-analytical behaviour can be seen to arise within Fermi 
liquid theory from the coupling between quasi-particles and incoherent spin fluctuations 
(quasi- holes in the triplet state), which give rise to non-analytical energy dependences in 
the self-energy 8J. These non-analyticities are in fact a very general feature of all normal 
Fermi liquids and their existence is not related to the particular details of 3 He systems . 

Hot symmetric nuclear matter (SNM), on the other hand, is an infinite fermionic system 
composed of nucleons at very high densities. In this ideal system, only the strong interaction 
among nucleons is considered and any other interaction, such as the electromagnetic one 
between protons, is neglected. The temperatures usually considered for this system are of 
the order of tenths of MeV and are somewhat small if compared with the typical energy 
scales of nuclear matter (the free Fermi energy at saturation density is about ep ~ 40 MeV, 
thus T jep ~ 0.5 for T = 20 MeV, the highest temperatures considered here). We expect 
that a finite temperature approach to SNM can fairly describe the properties of the hot 
environments that could exist either inside the cores of supernovae at the latest stage of 
their evolution [3] or in the collisions of heavy nuclei at intermediate energies [k|. The 
TD properties of these systems and, in particular, the entropy are important quantities 
for the understanding of astrophysical and heavy ion physics phenomena. In core-collapse 
supernovae, for instance, the evolution and dynamics occur at a fixed entropy per baryon 



ll| . Moreover, the entropy production in multi- fragmentation events in heavy ion collisions 



is considered to be a crucial quantity to determine the mass fragment distribution 121 ] . 

Hot SNM has been usually studied in a mean-field approximation with effective phe- 
nomenological nucleon-nucleon (NN) forces, such as the Skyrme or Gogny interactions 
[13I . [bjj ]. or within a Relativistic Mean Field approximation 15]. Other many-body ap- 
proaches that have been used in the study of SNM include lattice models 16) or three-loop 
calculations within chiral perturbation theory at finite temperature 171 ] . However, when 



dealing with realistic NN potentials, more sophisticated many-body techniques are needed. 
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The strong short range repulsion and the tensor components of these potentials modify sub- 
stantially the many-body wave function, which is not anymore well described in terms of 
a free Fermi gas Slater determinant. Particle-particle correlations, for instance, are crucial 
for a correct description of SNM properties from realistic NN potentials. The Brueckner- 
Hartree-Fock (BHF) approach accounts for such correlations by means of a summation of 
an infinite series of suitable ladder diagrams [18]. In fact, the BHF approximation arises 
from a well-defined expansion for the energy of a fermionic system at zero temperature, 
the so-called Brueckner-Bethe-Goldstone expansion [19(. At finite temperature, a similar 
summation can be achieved in the so-called Bloch-de Dominicis (BdD) approach 2(|. How- 
ever, this approach is no devised to reproduce the energy of the system. Instead, it aims 
to compute the TD grand-potential and, from it, all the TD properties of the system. In 
particular, the entropy includes the correlations embedded in the approach |2ll ] . 

Traditionally, however, the BHF approach has been extended to finite temperatures in a 
more naive way: the energy of the system is computed from a simple generalization of the 
T = formalism to finite temperature and the entropy of the system is computed from a 
mean-field expression 22j, |23[]. Finally, let us also note that relativistic BHF-type calculation 



24 



25|. 



at finite temperature have also been performed in the literature 

A consistent treatment of correlations in quantum many-body systems requires the inclu- 
sion of particle-particle and hole-hole scattering terms. The usual approach is the so-called 
ladder approximation or T-matrix approximation, in which an expansion of the single- 
particle Green's function in terms of diagrams is achieved. Such an expansion goes beyond 
the BHF approach in the sense that both the propagation of particles and holes are treated 
in the same footing. Within this approach, the propagator of a correlated fermion is not 
simply described in terms of one quasi-particle pole. Instead, the strength of a momentum 
state k is fragmented over a wide range of energies. The process of finding the Green's 
function for such a many-body system is unavoidably self-consistent, since the propagation 
of a nucleon is affected by the interactions with the surrounding nucleons, which in turn 
are also described in terms of Green's function. Although one can easily write down the 
diagrammatic expansion that gives rise to this Self-Consistent Green's Function (SCGF), it 
has taken a lot of time to have a numerical solution of the full problem. A major well known 
issue is related to the so-called pairing instability that appears in the zero temperature for- 



malism when the propagation of holes is considered 26|, |27|, [28] . Nevertheless, calculations 
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have been performed at the SCGF level at zero temperature within different approximations: 



with a quasi-particle self-consistent approach 



291 ]. with a discretized parameterization of the 



spectral function 30( or with a simplified separable NN interaction 31]. In all these ap- 
proaches the final output is the single-particle spectral function of a nucleon inside matter, 
which can also be obtained from other many-body approaches, such as the Correlated Basis 
Function theory 32]. 

In the case of SCGF, a finite temperature treatment has been the keystone for obtaining 
a complete numerical solution from a realistic NN potential 



33 



34 



35l. Once the self- 



consistent propagator is obtained in this method, one can easily obtain information on 
the microscopic properties (such as the momentum distributions, self-energies or spectral 
functions of the nucleon) as well as on the bulk properties of the system (the energy per 



particle via the Galitski-Migdal-Koltun (GMK) sum rule, for instance [36|, |37(). For a 
complete TD description of the system, however, one should compute the relevant TD 
potential of a statistical quantum mechanical system, i.e. the free energy. A suitable 
calculation of the entropy is thus required if this formalism is to be used in any practical 



in 



description of hot SNM. Here, we will follow the Luttinger-Ward (LW) approach [( 
which the grand-potential is computed from the full single-particle propagator. An analysis 
of the properties of SNM within this formalism has been recently published by Soma et al. 



39] . In the following we shall show that one does not need to compute TD quantities within 



a full LW approach as done in 39], provided that some approximations for the entropy are 
valid. Finally, let us notice that the LW formalism has been widely used in other many-body 
jhysics problems, ranging from relativistic plasmas {40I to resonances in heavy ion collisions 



4l|. 



In the following section, we will describe in detail the SCGF approach at finite tempera- 
ture and the LW formalism that we will use in our calculations of the entropy. The numerical 
results derived from this formalism will be divided in two different parts. Section III will be 
devoted to the microscopic results, while Section IV will describe the bulk TD properties of 
a correlated system of nucleons. Finally, a brief summary will be given in Section V. 
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II. FORMALISM 



A. The ladder approximation 

With the help of the single-particle propagator, we can obtain all the one-body (and even 
some two-body) properties of a many-body system 42|. In the Green's function approach, 
one aims to compute the single-particle propagator which, in the grand-canonical ensemble 
is defined according to: 

iG(kf,kY) = Tr {pT[a k (t)al(t')]} , (1) 
where we have introduced the density matrix operator: 

p = l-P&-tf) , (2) 

Zj 



and the partition function: 



Z = Tr\e-^-^Y (3) 



In Eq. (pQ), T is the time-ordering operator in such a way that the Heisenberg creation 
(annihilation) operator a k (t) (a k (t)) with the largest time argument t (or it in the case that 
t is imaginary) is put to the left, with a minus sign included for each commutation. In these 
equations (3 denotes the inverse temperature and /x, the chemical potential of the system. 
For simplicity, we will neglect the spin-isospin structure of this propagator in the following. 
Finally, the traces Tr are to be taken over all the energy and particle number eigenstates 
of the system. The cyclic invariance of these traces imply the following quasiperiodicity 
condition for the Green's function: 

G(kt = 0, k't') = -e^G{kt = -i/3, kV) . (4) 

One can thus Fourier transform the time dependence of the propagator in terms of some 
Fourier coefficients G(k, z v \ depending on the discrete Matsubara frequencies z v = ( - 2 ^- >7r + 
/i. A Lehmann decomposition of these coefficients is achieved by means of the spectral 
function A(k, lu): 
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The function G(k, z v ) evaluated at the Matsubara frequencies can be analytically continued 
for all non-real z. Using the Plemelj formula, the spectral function can be related with the 
values of G close to the real axis (we use the notation uj + = u + irj): 



A(k,u) 



-21mG{k,uj + ) . 



(6) 



The single-particle Green's function can be obtained from Dyson's equation, which for a 
translationally invariant system reduces to the algebraic equation: 

k 2 



u - E(k,u) 

2m 



G(k,u) 



(7) 



where Y,(k,u) denotes a complex self energy. By expanding the self-energy in terms of one- 
particle Green's function, it can be shown that £ and G share the same analytic properties. 
Thus, one writes the following spectral decomposition for S: 



^HF 



(*) + 



2tt z — uj 



where the width V(k,u) is related to the imaginary part of £ by: 



T(k,uj) = — 2ImS(/c, u 



+J > 



(8) 



(9) 



and it is real and positive for all energy and momenta [43(. The term Yi HF {k) in Eq. (jSJ) is 
a real energy-independent generalized Hartree-Fock contribution to the self-energy. 

The self-energy can be derived from an in-medium two-body interaction (the so-called 
scattering T-matrix) which includes the correlations induced by the strong short-range and 
tensor components of a realistic two-body NN force. Within the ladder approximation, one 



can indeed relate the self-energy and the retarded T-matrix by 44 



45j: 



/A 6 y r°° a,., i 
— J ^(kk'llm T{uo + uo' + )\W) A A{k\J) 

x [f{u')+b{u + u')\, (10) 



where we have introduced the Fermi-Dirac: 



_|_ I 



(11) 



and the Bose-Einstein: 



distributions. T he p ole of the Bose function b(Q) at Q = 2/x is exactly canceled by a zero 



in the T-matrix 



27 



461 ] and thus the integrand remains finite as long as the T-matrix does 



not acquire a pole at this energy. This pole appears only for temperatures below a certain 
critical temperature and it is closely related to the onset of pairing among nucleons 0, 4t| . 

Within the ladder approximation, the in-medium T-matrix is determined as a solution 
of the integral equation: 

(kk'|T(fi + )|pp% = (kk'|V|pp% 

x (qq'|T(ft + )|pp%, (13) 

where we have introduced the two-particle Green's function of two non-interacting but 
dressed nucleons: 

guk h, n + ) = r ^ r ^ A( kl , ^(h, ^-fw-fw . (14) 

Diagrammatically, the T-matrix approximation is depicted in Fig. [TJ where the dressed one- 
particle propagators are given by two double lines. In general, the equations of motion arising 
within the Green's function method couple the iV-particle propagator to the N + 1-particle 
propagator whenever we have a two-body interaction. The Dyson equation for the single- 
particle Green's function involves, for instance, the full two-particle Green's function Gu. In 
the Self- Consistent Green's Function approach within the ladder approximation, however, 
one ignores the effects of three- and more particle propagators by using a function, G° n , which 
is a product of two dressed one-body propagators. In that way, one solves self-consistently 
the equations for the one-body G and the approximated two-body G°jj propagators. 

If one wants to solve Eq. fflBl in an efficient way, a partial wave decomposition is needed. 
After expressing G° :I as a function of the total P = ki + k2 and the relative momenta 
k r = (ki — k2)/2, one can perform the usual angle-average approximation on the angle 
between P and k r , so that G° n is expressed only in terms of the modulus P and k r : 



g° h (p, k r , n+) = ^j l d(cos#) g° h (|p + k r | , |p - k r | , n + ) . 



(15) 



This approximation leads to a decoupling of the partial waves with different total angular 
momentum J, which in its turn implies that Eq. (TTOT) becomes a one-dimensional integral 
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equation: 

(k r \Ti!?(p,n + )\K) A = (h\vl[7\K) A 

x (^|T/fJ(^,n + )|^) A , (16) 
By summing over all partial waves, we get the T- matrix as needed in the solution of Eq. ( [131) : 
(kk'|ImT(fi + )|kk% = -i- Y (2J+1)(2T+1) 

47T z — ' 

J5TL 

x (g(k,k')|ImT/f T (P(k,k'),fi + )|g(k ! k'))A. (17) 
The only remaining piece is now the generalized Hartree-Fock contribution to the self-energy: 

^ HF ^ = h E( 2J+1 )( 2T+1 ) / ^(g(k,k')iF/f T ig(k,k')Un(A ; ') ) (18) 
where we have introduced the momentum distribution: 

<k)= r ^A(k,u)f(u). (19) 



oo 



2tt 



In a self-consistent procedure, Eqs. (IE]), (TTDTl . (jT3 T) - ( TT6T) are solved at a given temperature 
and density. The chemical potential /i is determined at each iteration by inverting: 

/d 3 k 
-j—^n(k,fi), (20) 

where v denotes the spin-isospin degeneracy of the system [y — 4 in the case symmetric 
nuclear matter). As a final output, once convergence is reached, one obtains the single- 
particle propagator from Dyson's equation and, from its imaginary part, the spectral function 
A(k, uj). From this function, several micro- and macroscopic properties of the system can be 
computed. For instance, the momentum distribution given by Eq. ( TT9"j) or the total energy 
per particle of the system, accessible from the GMK sum rule: 

where A is the total number of particles of the system. Further on, we will assess the 
question of how to compute the partition function and the entropy of a system of nucleons 
from the spectral function A. 
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Before doing so, however, let us consider some interesting approximations to the SCGF 
method. Using Dyson's equation, it is straightforward to show that the spectral function 
can be written as: 

Ak, = : mu) - r _„ , 12 • (22) 



2 



r(fc,t>)) 

2 



t is interesting to note that this is a positive defined function for all energy and momenta 



43j, fulfilling the following sum rule 48|, |49j, |50I |: 



/OO J 
^A(k,u) = 1. (23) 

One can obtain a simplified set of equations by taking the zero- width limit, r — > in the 
previous expression, thus obtaining the 5-function: 

A(k,u) = 2-k 5[u - e qp (k)] , (24) 

where the quasi-particle spectrum is derived by means of the real part of the self-energy: 

k 2 

e qp {k) = — + ReZ{k,e qp {k)). (25) 

This is the so-called Quasi-Particle Self-Consistent Green's Function method. In the fol- 
lowing, we will not use this approximation. However, we shall extensively use the 5-peak 
approximation to the spectral function, Eq. (1241) . with the single-particle energies e qp {k) 
defined by Eq. (1251) with the self-energy of our SCGF calculation. 

In addition, we will also consider the BHF approach to nuclear matter. This can be 
derived from the SCGF method if one takes the quasi-particle approximation, Eq. ( 1241 ). 
and, in addition, one makes the following substitution: 

(1 - /H - /(a/)) - [1 - f(u)\ [1 - /(a/)] (26) 

in the two particle propagator, Eq. (1141) . In the T — * limit, this expression becomes the 
Pauli operator and G° n reduces to the particle-particle propagator. From this approximation, 
one can see that the T-matrix reduces to the well-known G-matrix in-medium interaction: 

(kk'|G(ft + )| PP % = (kk'|y|pp% 

+ I A, /^(kk^lq^t 1 -^^^" 1 -^^^" 



(2tt)3 J (2tt)3 * Q + _ £BHF ( q) _ £BHF {q') 

x (qq'|G(fi + )|pp'>A. (27) 
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where the single-particle spectra of nucleons in the BHF approach, e BBF {k\ are given by 

k 2 

BBH F (k) = — + Re E BHF {k) . (28) 
2m 

In this approximation, the self-energy is given in terms of the G-matrix by: 

Z BH F(k) = j ^f[ £BHF (k')](kk'\G[Q = e BHF (k) + e BHF (k')]\kk') A . (29) 
and the total energy per particle is obtained from: 



E BHF y r d 3 k 



^-+ 1 -ReE BHF (k) 



A P J (2,)3 [ 2m ' 2 — ■^f[^HF(k)] (30) 

We shall make two more comments concerning this finite temperature generalization of the 
BHF approach. On the one hand, one must be aware that hole-hole propagation is the cause 
of the Bose term in Eq. (fIU|) . If the BHF is to be derived from the SCGF approach, one 
should neglect this term, so that the imaginary part of the self-energy reads: 

lmZ BHF (k) = y"^/[e B ^(A;0](kk , |ImG'[fi = ^F(A;)+^F(A ; ')]|kk%. (31) 

On the other hand, one should say that the BHF approach obtained from this approximation 
of the SCGF method is not fully justified from basic first principles. Indeed, the finite- 
temperature generalization of the T = Bethe-Goldstone expansion is given by the BdD 
approach. In that theory, a similar expression to Eq. (127D for the in-medium temperature- 
dependent interaction can be obtained, but in principle there is not a straightforward relation 
with Eq. fl30|) for the energy per particle of the finite temperature system. Nevertheless, one 
can see that this approximation is a reasonable one, because the dominant diagrams of the 
BdD expansion reduce to this finite temperature BHF approach at low temperatures 21]. 
In the BHF approach that we will use here, the free energy is obtained from the energy per 
particle of Eq. (I30I) together with the mean-field expression for the entropy of the system: 

S BHF = -uj -0^ [f[e B HF(k)] In f[e BHF (k)} + (1 - f[e BHF (k)}) ln(l - f[e BHP (k)])} . 

(32) 



B. Luttinger-Ward formalism 



From a statistical mechanics point of view, the macroscopic information of the system is 
contained solely in the partition function. Once this function is known, its derivatives give 
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access to the TD properties of the system. However, from a microscopic point of view, it 
is well known that the one-body properties can be derived from the dressed single-particle 
propagator of the system. Now one may ask whether there is a connection between both 
functions and, in particular, if the single-particle propagator is enough for building the 
partition function. The answer, given by Luttinger and Ward [38] more than forty years 
ago, is positive. The expression for the partition function thus obtained turns out to have 
some interesting properties that were later on exploited by Baym [5l| in his well-known 
discussion of the TD consistency of many-body approaches. The starting point is the so- 
called Luttinger-Ward expression of the partition function of the system: 

\nZ = TrV(k,z v )G(k,z J ,)+Trln[-G- 1 (k,z v )\-$[G], (33) 

where the functional has been introduced and where the trace Tr means a sum over 
all momentum states and Matsubara frequencies: 

Tr -> eZ " n > ( 34 ) 

with i] = + small and positive and such that liniR e r^Re z = oo. In TD equilibrium, the 
grand-partition function is stationary under variations of the Green's function: 

SlnZ „ , . 

= 0, (35) 

Go 



SG 

and thus the functional $ should satisfy the following condition: 



<5$ 
5G 



= £ . (36) 

Go 



To obtain the previous expression, we have used Dyson's equation, Eq. (J7J) in order to 
functionally derive the self-energy £ with respect to G. In fact, within the LW formalism, 
£ is a functional of the propagator and one can take Eq. (1361) as its definition. Usually, 
however, in a given many-body approach one assumes a certain approximation to the self- 
energy (in our case, the ladder approximation). The functional $ is then fixed by Eq. fl36|) 
and the propagator is given by Dyson's equation, G^ 1 = Gq 1 + E. These ideas can be 
depicted diagrammatically. For the ladder approximation, we show in Fig. [2] the n-th order 
contribution to the self-energy £( n ) (where n denotes the number of bare interaction lines 
in the diagram). By closing this diagram in its free vertices with a propagator, we obtain 
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the corresponding n-th order contribution to $. At each order, however, this contribution 
should be divided by a factor In that takes into account the 2n possible places where one 
can cut each propagator to obtain E from Eq. (155]) [51]. The corresponding series for the 
functional is shown in terms of diagrams in Fig. [3J 

The evaluation of the sum in the partition function has to be done with special care 
because of the cut in the logarithm. The final expression for the grand-potential Q = —T In Z 
reads: 

/d 3 k f°° duj ( ^ 

— J —fXu)2l^\n[-G-\k,u + )]+nk,uJ + )G(k,u + )^ 

+ T$[G] . (37) 

In this expression, the propagator and the self-energy are computed above but close to the 
real axis. We have also introduced the Fermi-Dirac distribution f(u), Eq. (TTTj) . One can 
now readily obtain the entropy by means of the TD relation: 



_ on 
b ~ df 



(38) 



The stationarity of Q with respect to changes in G is now very useful, because it implies 
that in Eq. ( 1371) only the temperature derivatives of the Fermi functions are needed: 



S = V SW? [^ dJ ^^[\n[-G'\k,^)\+T,(k,^)G(k, 



+ §f T ^ G \ ■ (39) 

This expression gives the entropy per unit volume of a correlated system of fermions as a 
function of G, S and $ and it is the fundamental equation from which we will derive most 
of our results. The usefulness and applications of this formula in the context of Fermi- 
liquids were extensively discussed in the pioneering work of Carneiro and Pethick |f|. In 
the following, we will closely follow this reference and discuss some of the more relevant 
approximations for our case. 

One can explicitly compute the imaginary parts of the two terms inside the integrals. 
The first term is the imaginary part of a logarithm, i.e. the phase of its argument. We can 
thus decompose the argument in its real and imaginary parts: 

k 2 i 

G-\k, uo+) = uo Re E(k, u) - ihnE(k, u+) = Re G~\k, u) + -T(k, u) . (40) 

2m 2 
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We consider that the logarithm has a cut in the real negative axis and we work in the 
sheet where ln(l) = 0. The arctan(z) function goes from — | to | and is normally used to 
obtain the phase of a complex number. Nevertheless, whenever the real part of G~ x becomes 
positive and provided that F(k,u) is positive defined, the complex phase of the logarithm's 
argument will be in the III quadrant of the complex plane and thus a factor 7r needs to be 
subtracted to the the arctan function in order to match the argument of —G^ 1 . The result 
of the imaginary part can thus be casted in the following form: 



Imjln [- G-\k,u + )]^ = arctan X(k, uj) - ttO [ReG~ 1 (k, < 



co)] , (41) 



where we have introduced the function X(k,u): 

which is nothing but the quotient between the imaginary and the real part of the complex 
inverse propagator. On the other hand, the imaginary part of the second term in the integral 
of Eq. ( 1391) is given by: 



S^{k,LU + )G(k,U + )j 



Im \ uj+)G(k, uj+) \ = ReYi(k,u} + )lmG(k,u} + ) + ImS(fc, u + )ReG(k,u + ) 

= -~ReZ(k,u)A(k,uj) - ~F(k,co)ReG(k,uj) . (43) 
Using Eqs. (HI!) and (j4"3"j) . we can divide the entropy S in two terms: 

S = S DQ + S' , (44) 



given by: 



and 



* - ->^Sw?L^ A{kM ^ (46) 

where we have introduced the function S(k,uj): 

E(k,u) = 27i9[ReG-\k,uj)] - 2 arctan A(A;, cu) + T(k,to)ReG(k,uj) . (47) 

S D ® is a dynamical quasi-particle (DQ) entropy which partially takes into account the cor- 
relations of the dressed particles in the medium. It includes some finite width effects, as 
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seen by the fact that it is computed with a non-zero T. The term S', on the other hand, 
accounts for higher order correlations. As it was shown in |6j, this term arises from the 
cancellation between the temperature derivative of T$ and the second term of Eq. (1461) . 
The only non-zero contributions that survive this cancellation come from terms in the per- 
turbation expansion that have at least two vanishing energy denominators. In the following, 
we shall make the assumption that S' is negligible. In this way, our formalism is simplified 
because there is no need to evaluate the $ functional. This assumption, of course, needs to 
be validated, and that is what we will do in the final part of this work. As we will see, our 
approximation leads to TD consistent results, which confirms that the contribution of 5" is 
small in the density and temperature range explored as far as short range correlations are 
concerned. In addition, it is worth to mention that it is precisely from the S' contribution 
that all the anomalous temperature dependences of the entropy arise. These anomalies, 
however, are mainly generated by terms in S' involving long-range correlations, which we do 
not consider in our approach. Thus, by restricting ourselves to computing the entropy with 
S D ®, we will loose these contributions and only find analytical (S ~ T, T 3 ) temperature 
dependences. 

To study the DQ entropy, it is interesting to analyze the properties of the function 
E(k,uj). The first term of H in Eq. (T4T1) is a step function with argument ReG -1 . For a 
fixed momentum k and as a function of the energy, this argument is negative for energies 
below the quasi-particle peak e qp (k), whereas it is positive for u greater than e qp {k). Thus 
the step function can be rewritten: 



At a fixed momentum, then, Hi equals zero at energies below the quasi-particle pole and 2ir 
above it. By using the relation: 



uS) = 2tt6 uj 



2m 



Re to) 



2n9 [< 



u - e qp (k)] . 



(48) 





(49) 



dT 



duo 



where we have introduced the function: 



a 



(«) = - f(u) ln/M + [1 - /H] In [1 - /(« 



)] 



(50) 
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the contribution of Si to the entropy is given by: 

Sir 1 = v \ - — — / 2-K0 \u - e m (k)\ = v \ - — - / duo— — = 

1 J (2tt)3 2tt du 1 qpK Ji J (2vr)3 du 

iVk a[e qp (k)]=S QP - (51) 



(2tt) 3 

This expression corresponds to the entropy of a system of undamped quasi-particles with 
real quasi-particle energies given by Eq. (T25]) . Whenever quasi-particles have long lifetimes, 
we expect it to be a good approximation to the entropy. Indeed, for any many-body ap- 
proximation where the quasi-particle energies are real (such as the Hartree-Fock case, for 
instance) the full DQ entropy is simply given by Eq. (|5ip . The rest of the terms in S can 
be rewritten as a function of \(k,u): 

E 2 (k,u) = — 2 arctan [\(k, lo)] , (52) 

for the second term and: 

for the third one. Their total contribution to the entropy is then given by: 

This expression involves a non-vanishing width T and it can thus be thought as a lifetime 
correction to the DQ entropy. It is clear that for infinitely long lived quasi-particles (r=0), 
this contribution will be zero, but for large widths it can have a non-negligible effect on the 
total entropy. 

Let us go back to Eq. fH5]) for the DQ entropy. After a partial integration and using 
relation we see that the following expression for the DQ entropy holds: 

sD0 -/S£^» B( ^»- (55) 

provided that the B spectral function is defined as: 

dZ(k,uj) 

B{k, uj) = — — — . (56) 

The expression Eq. ( |55|) has several interesting properties. First of all, in the free and the 
Hartree-Fock cases the B function reduces to a delta peak and S D ® becomes the expected 
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expression for the entropy per particle. However, this does not mean that S D ® neglects the 
finite width of quasi-particles, as we already commented. In addition, it is easy to check 
that the B function fulfils the following sum rule: 



f°° dec „„ , 1 , 



1 oo 
I — oo 



1 . 



(57) 



Finally, let us remark that this expression is somehow intuitive, in the sense that it is a 
product of the statistical weighting factor for the entropy <j(uj) times a spectral function 
that takes into account the width of quasi-particles. 

An alternative way to obtain the entropy within the LW formalism is obtained by starting 
from Eq. ([57)1 and using the fact that: 







/3f(u) = -— In [1 + 
ouj - 



Integrating by parts one easily obtains: 

d 3 k 



n 



where we have introduced the B function: 



B(k,u) = -^-2Im|ln [- G-\k,u + )] + £ 



(k,LU+)G(k,uj + ] 



(58) 



(59) 



(60) 



The entropy is now obtained via the temperature derivative of Eq. (I59p which will only affect 
the explicit temperature factors thanks to the stationarity condition Eq. ( 1351) : 



S 



v 



d 3 k 



00 A, ; f) 

a(u)B(k,u) + -^-T<S>[G}. 



2tt 



dT' 



(61) 



Different expressions for B can be obtained depending on whether the derivative or the 
imaginary part are taken first. By taking first the imaginary part and then the derivative, 
we see that the following relation holds between 5 and B: 

B(k,u) = —E(k,u) + ^ ' J A(k,u) + ReE(k,oj) . 62 

OUJ OUJ OUJ 

For the first term of Eq. ( 160]) . however, it is instructive to follow the opposite direction, 
deriving first and taking after the imaginary part we obtain: 



B x %u>) = -21m |^-ln [ - G~ 1 (k,uj + )] j = -21m | 



A(k,u) 



1 - 



dReZ(k,uj) 
duj 



Re G(k, uj 



2lm tG(k,w+) 
dT(k,u) 



d"E(k, uj^ 
duj 



du 



(63) 
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The derivative of the second term is easily computed and can be separated in two parts: 



B 2 (k,u) 



dReG(k,u) dT(k,u) 
-T(k,u) + Re&(k,u) 



du 



du 



and 



B 3 {k,u) 



dReT,(k,u) ,., , .dA(k,u) 
-A(k,u)+ReE(k,u)- 



du 



did 



Now we can write a compact expression for the B function: 
B(k,u) 



... , dA(k,u) n . dReG(k,u)^ /7 . 
A(k,u) H K -^—^ReZ(k,u) H v ; F(k,u) 



(64) 



(65) 



(66) 



du du 

In addition, if one uses the fact that B% equals the last two terms of Eq. (1621) . it is easy to 
obtain the following relation between the B and the B spectral functions: 



B(k,u) = B 1 (k,u) + B 2 (k,u) = A(k,u) 



dReZ(k,uj) 
du 



dReG(k,u)^., , , , 
+ ^^T(k,u),(Q7) 



which also gives B as a function of A and the real and imaginary parts of the propagator. 

The expression for the entropy of Eq. flBTT) . in terms of the $ functional plus a term 
containing the integral of a statistical factor u{u) and the weighting function B(k,u) has 



already been obtained in, for instance, 39(. However, one should take into account that 
our approximation for the entropy, S D ®, differs from the first term of Eq. (I6ip . In our case, 
we have neglected the terms of Eq. (|46p which are in fact canceling each other to a certain 
degree, whereas by approximating the entropy with the term of the B spectral function 
one would be ignoring this cancellation. In fact, there is no reason to believe that the first 
term of Eq. flBTj) should be a good approximation to the full entropy, while the DQ entropy, 
given by the convolution of the statistical factor o{u) and the B spectral function, gives very 
reasonable results as we shall see in the following. 

Also illustrative is the following decomposition of S D ® into two terms: 

d 3 k r°° du 



(27T) 



2tt 



cr(u) A(k, u>) 



and 



v 



d 3 A; 
(2^)3 



du; , . f dReG(k,u)^, 7 , .., .dReT,(k,u) 
a M \ a F{k,u) - A{k,u) 



2tt 



du 



du 



(68) 



(69) 



This justifies somehow the naive generalization to the expression of the entropy that has been 



used in the literature 



52 



53| and which consists in approximating the entropy by formula 



17 



(!68|) . This is, of course, not justified from TD grounds, but it would be a natural extension 
of Eq. ffl9l) to the case of the entropy. In particular, within a quasi-particle approximation, 
when both A(k, u>) and B(k, uj) become quasi-particle delta functions, S DQ and S± coincide 
and become the quasi-particle approximation to the entropy, S® p , 

/r\ 3 k 
(2^)8 ^M*)l- ( 7 °) 

The DQ entropy, however, goes beyond the naive quasi-particle approach. It introduces the 
corrections of Eq. ( 1691) which, as we will see later on, are non-negligible. 

Another justification for Eq. ( |55i) comes from the generalization of the ^-functional tech- 



nique to non-equilibrium quantum systems. In Ref. [54j, it was shown that, within certain $- 
derivable approaches out of equilibrium, an if -theorem could be proved for a non-equilibrium 
kinetic entropy expressed in terms of the full Green's function and the self-energy. When 
dealing with equilibrium systems, this kinetic entropy reduces to the sum of S D ®, the local 
or Markovian part of the kinetic entropy (Si oc in the language of Ref. (541]), plus S', the 
memory or non-Markovian part of the entropy (S mem ), which coincides with the expression 
of the entropy Eq. ( 144)) . 

If S D ® is close to the real entropy of the correlated system S, we should recover TD con- 
sistency, i.e., the microscopically computed quantities should coincide with the macroscop- 
ically com put ed ones. The ladder approximation is known to be a well-defined ^-derivable 



approach 



5l| . In particular, the (microscopic) chemical potential fi computed from the 



normalization condition: 

f d 3 k f°° du , „, 

should coincide with the (macroscopic) chemical potential coming from the TD expression: 

dF(p,T) 



/i(p,T) 



(72) 



dp 

with F the free-energy per unit volume. At T = 0, this has been numerically checked for 
the SCGF with separable potentials in [55(|. However, it is also well known that some many- 
body approximations do not fulfil this check of consistency. The BHF approach, for instance, 



badly violates the Hugenholtz-van Hove theorem [56( which states that, at saturation density 
the chemical potential pL and the free energy per particle FJA should coincide. The difference 
between these two quantities can be as large as 20 MeV [571 ] . 
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In the following, we will present the numerical results for the B spectral function and 
the different approximations to the entropy obtained from SCGF's calculations. Our aim is 
twofold. On the one hand, we will show that the computation of the DQ entropy is enough to 
maintain TD consistency; that is, that for the range of densities and temperatures considered 
gDQ gj ves a f ree energy that respects the Hugenholtz-van Hove theorem. This consistency 
is embedded in the ladder approximation (as shown by Baym [5l|), but it is lacking in other 
many-body approaches. In addition, the TD consistency of our results seems to indicate 
that the S' contribution to the entropy is not crucial in the nuclear matter case. Note that 
an exact evaluation of S' would need the knowledge of the ^-functional. This functional has 

n 

recently been computed also within a ladder approximation [39J, in a SCGF computation 
that differs from ours only in some numerical details. The results that we will present 
(specially those concerning the temperature dependence of the entropy) agree substantially 
with those of Ref. |39j ]. We believe that this is an indication of the consistency of both 
approaches. 



III. MICROSCOPIC RESULTS 



All the results quoted in this and the following sections have been obtained with the finite 



temperature SCGF approach of Ref. 33| using the CDBONN potential £ 



. In the numerical 



treatment, partial waves up to J = 8 have been included. The Born approximation has been 
used for J > 3. The quoted BHF results have been computed with the same NN potential 
with partial waves up to J = 4. None of the calculations includes three-body forces. We 
are thus not able to reproduce the saturation point of SNM. In this sense, the results here 
presented should be taken as a first study of the TD properties within SCGF theory, focused 
on the effects that correlations induce on the entropy of SNM. 

In the previous section we have mentioned that the properties of the B spectral function 
are very close to those of the usual spectral function A(k,u). It fulfils a sum rule and it 
accounts somehow for the effect of the width of quasi-particles in the DQ entropy. It is thus 
natural to compare the two functions in the same plot. We can get a rough idea of the 
differences of both functions following an argument first proposed by Carneiro and Pethick. 
Let us express the spectral function A as a function of the real and imaginary parts of the 
self-energy. We will of course obtain the well-known Lorentzian-like function of Eq. (122ft . 
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For a given momentum, the spectral function will have a peak around the quasi-particle 
energy of height A ~ 4/T(fc, e qp {k)). On the other hand, the B spectral function can also be 
rewritten in terms of the self-energy. Starting from Eq. (14 7p and taking the derivative with 

respect to the energy, we get: 

1 T 3 {k,u) 



B(k,u) 



[w-g-IteE&a; 

r 2 {k,u 



r(fc,ai) 
2 



1 2 



-i 2 



[w-^ -Re £(*,«;)] + 



r(fc,ai) 



1 



OJ 



dReZ(k,u) 
du 



k 2 AdT(k,u), 
— -ReZ(k,u)j — ^(.73) 



<9c<j 



If one assumes that the frequency dependence of T and Re £ are smooth close to the quasi 
particle energy, we will have: 

i r 3 {k,u) 



B(k,u) 



"r(fc,a;)" 


2" 


to 





2 ' 



(74) 



which corresponds to a function which decays faster than a Lorentzian close to e qp (k), but 
which has a stronger peak at the quasi-particle energy, B ~ 8/T(k, e qp (k)). 

One can check that this schematic scenario is true in Fig. HJ where we show the B (full 
lines) and the A (dashed lines) spectral functions as a function of the energy at three different 
momenta k = 0, k = kp and k = 2kp, at the experimental saturation density p = 0.16 
fm -3 and at a temperature of T = 10 MeV. In all the three panels, corresponding to the 
three different momenta, we see that both functions are peaked around the same energy 
values, corresponding to the quasi-particle energies given by Eq. (|25|) . The peaks shift from 
negative values (with respect to the chemical potential) to positive values when going from 
zero-momentum to higher momentum states, just following the position of the quasi-particle 
peak. However, while the A spectral function has high-energy tails that contribute in a non- 
negligible way to the total strength of the nucleon, the B spectral function has lower and less 
extended energy tails. This is easily understood if one considers that both functions fulfil 
the same sum rule. Since the B function has a higher quasi-particle peak, the strength of 
the peak is contributing substantially to the total sum rule and there is no need to generate 
high-energy tails. The presence of these high-energy tails in the A function is an indication 
of the importance of the correlations that go beyond the mean-field approach [331 ] . Thus, 
the lack of such tails in the B function is signaling somehow that these correlations will have 
a small influence in the total entropy of the system. 
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This idea is also in accordance with the behaviour of the width of both spectral functions. 
Far away from the Fermi momenta, both functions are relatively broad around the peak. 
Again, in the case of the A function this is a consequence of the correlations that redistribute 
the nucleon single-particle strength within a wide range of energies. The B function has a 
smaller width, which indicates that it is less affected by correlations. Close to the Fermi 
momentum, however, both functions approach a delta-peak behaviour, reminiscent of the 
fact that at zero temperature, even when correlations are included, the spectral function 
has a delta-peak contribution. At this momentum and for the temperature considered, the 
B function is narrower and much more peaked than the usual spectral function A. 

Also note that the values for the B function are positive for all the energies and momenta 
that we have considered in our investigation. This is in contrast to the weighting function 
B, which is defined in Eq. ( 160]) and has been used in Ref. [391 ] . The fact, that the evaluation 
of the entropy using the weighting function B exhibits strong cancellation effects (see Fig. 4 



of Ref. 



39(| ) may be taken as an indication that the splitting of the entropy into the two 



contributions according to Eq. ( ISTT) might not be optimal. 

In order to understand the density dependence of the DQ entropy, we show in Fig. [5] the 
B spectral function as a function of the energy for different densities (p = 0.1, 0.2, 0.3, 0.4, 0.5 
fm -3 ) at the same three different momenta previously compared and at a fixed temperature 
of T = 10 MeV. In addition, we plot with a dotted line the statistical weighting function a(uj). 
It is precisely the product of these two functions, integrated over energies and momenta, 
that gives rise to the DQ entropy, so it is interesting to study their overlap. 

The general features of the B spectral function as a function of density are very close to 
those of the usual spectral function A. In the case of k = 0, the quasi-particle peak moves to 
more and more attractive energies as density increases, reflecting the fact that the binding 
energy of a zero-momentum nucleon increases with density. Above the Fermi surface (at 
k = 2fcp), the situation is the opposite and the peak of the B function moves to higher 
energies with increasing density. The width of these peaks, both at zero momentum and 
at twice the Fermi momentum, is broadened with density. This is in accordance with the 
naive idea that correlation effects increase with density. In addition, as a consequence of 
this broadening, the high energy tails (visible at high positive energies for the k = state) 
decrease with density, allowing the sum rule Eq. floTl) to be fulfilled. 

The situation is different at the Fermi surface: when the density is increased, the peak 
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remains at a fixed energy oj = //, while its width becomes narrower and concentrates more 
strength. This can be understood if one takes into account that, as already commented, at 
zero temperature, i.e. for the fully degenerate system, the correlated A spectral function 
shows a delta peak which would also be present for the B spectral function. At a fixed non- 
zero temperature, however, the system moves towards the degenerate limit (the ratio T/ep 
decreases) with increasing density and thus the B spectral function becomes closer to a delta 
peak. This is actually what can be seen in the central panel of Fig. [51 At high densities 
(p > 0.2 fm~ 3 ), a clear separation between the quasi-particle peak and the background 
contribution to the B spectral function is observed. 

It is clear from Fig. [5] that the quasi-particle peak and the peak of the a function only 
coincide for momenta close to kp and energies around u = \i. Thus, the important contribu- 
tions to the DQ entropy of the system will be that of the momenta close to the Fermi surface 
and the energies close to the chemical potential. It is precisely the interplay between a and B 
that gives rise to the density dependence of the entropy. Since the value of B at k = kp and 
u — fi increases with density, one may expect that the entropy per particle would increase 
with density. However, it is also true that, for lower densities, the quasi-particle peak is 
closer to fi at all momenta and thus there are contributions of the quasi-particle peak for 
momenta not necessarily close to kp. In fact, when these contributions are summed up, one 
finds that the entropy per particle decreases with density. 

To illustrate these results, we plot in Fig. [6J the momentum-dependent integrand of 
Eq. (ESJ): 

/CO 1 
=-a(u)B(k,u), (75) 
-oo 

which measures the contribution to the DQ entropy from each momentum state. It is 
clear that, as density increases, the integrand becomes larger at the Fermi surface but less 
extended in momenta. This is in agreement with the previously discussed ideas, i.e., that for 
less degenerate systems the contributions at all momenta are relevant, while for degenerate 
systems the contribution of the k = kp state is the most important one. In addition, the 
circles in the figure show the ( function obtained within a quasi-particle approximation: 

C QP (k) = ^k 2 a[s qp {k)\ , (76) 

with the quasi-particle peak given by the maximum of the spectral function, Eq. (125]) . The 
differences are only relevant for the lowest densities and in a range of momenta close to the 
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Fermi momentum. This is again a signature of the small role played by the correlations that 
deplete the quasi-particle states on the entropy. Therefore, we expect that the quasi-particle 
approximation to the entropy, Eq. (1701) . will describe the full DQ entropy S D ® quite well. 

In order to gain insight in the temperature dependence of the DQ entropy, we show in 
Fig. [7] the B spectral function as a function of energy for the same three different momenta 
considered previously at a fixed density p = 0.16 fm~ 3 and at 5 different temperatures 
T = 4, 8, 12, 16, 20 MeV. It is clear that, for all momenta, the variations of temperature 
mainly result in changes of the width of the quasi-particle peak, while the position in energy 
of this peak relative to the chemical potential is almost not changed. In addition, the 
momentum states far above the Fermi surface are not affected by temperature, as it is seen 
in the lowest panel, corresponding to k — 2kp. At the Fermi surface, on the other hand, the 
effects are much more important. As temperature increases, the height of the quasi-particle 
peak decreases, while its width increases. Moreover, at the lowest temperature (T = 4 
MeV) a clean separation is observed between a quasi-particle peak and a particle and hole 
background. This separation is softened at T = 8 MeV and disappears completely above 
this temperature. Such a behaviour is again understood in terms of the degeneracy of the 
system. The lower the temperature, the higher the degeneracy and the smaller the width of 
the B spectral function. For the k = state, a similar situation is found. The peak lies below 
the chemical potential, and it is clearly splitted from the particle background at T = 4 MeV. 
For temperatures above T = 12 MeV this separation disappears and a smooth transition 
from particle to hole states is found in the B function. It is also interesting to notice that the 
width of the peak remains more or less constant, thus indicating that temperature- induced 
effects on the width of the quasi-particle peak are more important at k = kp. 

As for the total contribution to the DQ entropy, the convolution between o{iS) and B(k, uj) 
is again crucial. At low temperatures, <j(uj) is very peaked around u ~ fi. The convolution 
will thus only be different from zero whenever the quasi-particle peak is close to fi, i.e., 
at k ~ kp. On the other hand, at higher temperatures ct(uj) is different from zero in a 
wider region of energies, which results in a non-zero convolution at all momenta. When we 
integrate over momenta, the final DQ entropy is higher for the higher temperature. Thus, in 
accordance with intuition, the entropy of this correlated system will grow with temperature. 
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IV. MACROSCOPIC RESULTS 



In this section we will explore the density and temperature dependences of the entropy 
computed within different approximations. To begin with, we show in Fig. [H] the density 
dependence of the entropy per particle at a fixed temperature T = 10 MeV. The approxi- 
mations to the entropy per particle that appear in this Figure are: 

• S DQ , the full dynamical quasi-particle entropy of Eq. (j4"5]l and equivalently Eq. (I5"5"|) 
(solid lines). 

• S QP , the quasi-particle approximation to the entropy of Eq. (1701 (dotted lines). 

• S BHF , the Brueckner-Hartree-Fock entropy computed from expression Eq. fl32|) 
(dashed lines). 

• S J 4 , the contribution to the entropy due to the A spectral function (dash-dotted line). 

As a general feature, we can say that all of these entropies decrease substantially with 
density, from values of around 2.5 at densities of around 0.02 fm~ 3 down to values of around 
0.4 for the highest density here considered, p = 0.5 fm~ 3 . 

One important result that arises from Fig. [8] is the fact that, at all densities, S D ® and 
S® p are very close. This is somehow in agreement with the idea that the inclusion of the 
width of quasi-particles has a small effect in the entropy. As discussed in relation with Fig. [HI 
the effect is larger at lower densities, where both approximations differ more, but it is never 
higher than a 5%. At high densities, the difference is so small that it cannot be appreciated 
in the Figure. This result is not at all intuitive. It indicates that S^, which is nothing 
but the difference between S D ® and S QP , decreases with density. But, since we have argued 
that S2 represents somehow the finite width effects on the entropy and since correlations 
grow with density, we would also expect it to grow with density. However, we have also 
seen that the higher the density, the smaller the width of the B spectral function at k = kp 
(which is the more relevant contribution at high densities) and thus the lower the effects of 
correlations. This is why at higher densities both approximations to the entropy tend to be 
similar. The fact that S®® is negative at intermediate densities (say from p = 0.05 fm~ 3 to 
0.30 fm~ 3 ) is quite interesting: in addition to stressing the fact that finite lifetime effects 
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to the entropy are small, one can say that it looks like correlations (i.e. the width of the 
quasi-particle) tend to order the system. 

The effects of the self-consistent propagation of holes are responsible for the difference 
between S D ® and S BHF . This can be clearly seen by rewriting Eq. (13"2"1) to give: 

S BHF = v J ^a[E BHF {k)\, (77) 

which differs from Eq. (ITUj) in the position of the single-particle peaks (given by s qp (k) in 
the first case and by EBHF(k) in the second) and also in the different values of the chemical 
potentials in the statistical factor a. Since the effect of the width on S D ® is small at this 
temperature, the difference between both entropies arises from the different quasi-particle 
energies and chemical potentials of the two approaches. Each of these quantities can differ 



by at most 20 MeV 



331 ] . However, the correction in the entropy is small probably due 



to a cancellation between both differences in the argument of a, where e(k) and \x are 
subtracted. In the intermediate density region, the BHF entropy has values that are about 
10 % below the DQ one. The presence of hole-hole correlation, thus, increases the entropy, 
i.e. the thermal disorder. This is of course related to the fact that hole-hole correlations 
tend to increase the density of single-particle states close to the Fermi energy. If one tries 
to parameterize the quasi-particle spectrum close to /i in terms of an effective mass m*, one 



33|. 



obtains larger values for the parameterization of the SCGF spectrum than for BHF 

Finally, we also show the contribution of S± to the DQ entropy. As we have already 
mentioned, this expression comes from a naive generalization to incorporate width effects, 
but nevertheless it gives a reasonable first guess to the entropy per particle. Intuitively, 
one would expect that, since the A spectral function is wider than the B one, the overlap 
between A(k, u) and <t(uj) at a given momentum should be higher and the final entropy 
would overestimate S D ®. However, this is not the case, except for the lowest densities. 
This can be understood from the the height of the quasi-particle peak for A being, roughly 
speaking, a factor of 2 lower than that for B. Thus, although more extended in momentum, 
the ((k) function for the A spectral function is smaller and gives rise to a lower entropy. 
The difference of both entropies is between 20 — 30% for the intermediate density region. 
The origin of such differences is the contribution of Eq. (1691) . which is the integral of 
two terms. Both terms are of the same order at p ~ 0.1 fm -3 but, while the contribution 
proportional to A decreases with density, the one proportional to T increases, and above 
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p = 0.3 fm~ 3 it carries more than 80% of the total correction. 

In Fig. [9] we show the temperature dependence of the entropy for a density p = 0.16 
fm -3 computed with the same approximations as discussed above. In addition, we have also 
computed another approximation, S NK , which is displayed with a double-dotted dashed 
line. This corresponds to the entropy of the correlated momentum distribution, Eq. (1191) . 
obtained from the mean-field- like expression: 

S NK = -v J -0^ |n(fc) In n(k) + [l - n{k)\ In [l - n{k)\ J . (78) 

One can say that all the approximations to the entropy of Fig. [9] decrease with temperature 
(as expected) and approach a linear dependence at low T. It is a well known feature of 
Fermi liquids that the slope coefficient for such a linear behaviour is proportional to the zero 
temperature density of states computed at the Fermi surface, N(T = 0): 

% = ^iV(0)T. (79) 

Each of our approximations goes to the T = limit with different slopes, and we can thus 
obtain different densities of states. To calculate N(0) for each approximation, however, we 
should extrapolate our results to the T = limit, which we cannot do within our approach 
reliably. Instead, from a low temperature expansion of expressions Eqs. fl55|) and fl68l) we 
obtain an expression of iV(0) in terms of the B and A spectral functions that we extend 
to finite temperatures. Namely, we get the "density of states" related to the B spectral 
function: 

/d 3 k 
——B(k,uj = p), (80) 

and the usual one, related to the A spectral function: 

N A (T) = u f -^A{k,u = ri, (81) 

where T denotes the fact that these have been computed at the finite temperature at which 
the spectral functions have been computed. Note that at low T and u = p the functions 
A and B differ basically by a factor Z (where Z denotes the renormalization of the quasi- 
particle pole): 

A(k, uj = n) ~ Z^(k, uj = p)B(k, lu = p) (82) 
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and thus the two densities of states will also differ by such a factor. Nevertheless, the 
density of states that gives the correct linear fit to the DQ entropy is that of Eq. (15U1) . In 
fact, we have numerically checked that it is this quantity that reduces to the well-known 
quasi-particle expression: 

Nqp{0) = (83 ) 
at low enough temperatures, where the effective mass is obtained through the derivative: 



m 



1 / de qp (k) 



(84) 



m 2m \ dk 2 

evaluated at the Fermi surface using the finite temperature SCGF quasi-particle spectrum, 
e qp (k). This can be seen in TableHJ where we give the densities of states computed with the A 
and the B spectral functions together with that obtained from the QP expression, Eq. ( 1551) . 
at p = 0.16 fm~ 3 for low temperatures. The effective is given in the fifth column 

of Table [H The numerical values confirm that at low temperatures the density of states 
from the B spectral function reduces to the quasi-particle one, hence indicating that from a 
TD point of view this is the correct density of states. However, from a microscopic point of 



view, the A density of states is the one which has been commonly used [59|, |60|]. One should 
keep in mind that in a mean-field approximation both of them would reduce to the same 
expression, Eq. (1551) . It is also interesting to note that if we compute the entropy with the 
help of the "densities of states" coming from the B spectral function, Eq. (1501) . and from the 
QP expression, Eq. (1551) (both of them evaluated at the corresponding finite temperature), 
and we use them instead of N(Q) in Eq. ( |79l) . we can reproduce the DQ expression of the 
entropy for temperatures up to T = 10 MeV with less than a 10 % discrepancy. 

In the context of nucleus-nucleus collisions at intermediate energies, there exists a grow- 



ing amount of experimental data 



61 



62] which should be useful to constrain the thermal 
properties of nuclei and nuclear matter. In particular, the liquid-gas phase transition and 
the caloric curve give a hint on the properties of nuclei at low temperatures. In the study of 
the caloric curve, it is customary to parameterize the excitation energy at low temperatures 



in terms of the so-called inverse level density parameter K (see [63(| for a theoretical descrip- 
tion), which is inversely proportional to the density of states introduced here. The values we 
obtain for K (defined as K~ l = ^N B (0)) are close to the Fermi gas value K ~ 14.6 MeV for 
p = 0.16 fm -3 . This can be understood from the fact that N B reduces to the quasi-particle 
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value of Eq. (IBB1 which, in addition, is similar to the free Fermi gas value because in our case 
the effective mass is almost equal to the bare nucleon mass at low temperatures. A word of 
caution must be raised, however. Our value for K is obtained from a calculation in infinite 
isospin symmetric matter in which only short range correlations are treated. Nevertheless, 
it is clear that a study of the inverse level parameter should include both the effects of finite 
size and long-range correlations, which are very important in determining the low energy 
excitations of nuclei. 

In the following, we compare the different approximations to the entropy that appear in 
Fig. [9j As a general trend, we observe the same features that we have already discussed when 
we have commented the density dependence of the various approximations to the entropy. 
The quasi-particle approximation S® p using SCGF energies reproduces the DQ entropies 
at all temperatures very well, especially below T = 10 MeV. The finite temperature BHF 
entropy describes the entropy of the system, with an error of about 15 %. This difference is 
thus quite small, which is again a signal that both the depletion of single-particle strength 
and the exact position of the quasi-particle peak are not that crucial in the final result of 
the entropy per particle. Let us also note that the inclusion of the hole-hole propagation in 
the quasi-particle peaks tends to increase the entropy. 

The A spectral function contribution to the entropy is again correct in the temperature 
behaviour, but still it gives a too small value for the entropy per particle of the system, 
with errors as large as a 30 %. As we have previously discussed, this is due to the lower 
quasiparticle peak of the A spectral function, which makes the S± contribution to the 
entropy lower than the S D ® entropy. The difference between these two functions is given by 
the S2 contribution to the entropy. This is composed of the two terms in Eq. ( |69|) . which 
have different relative weights as temperature changes. While the term proportional to the 
A spectral function amounts to 90 % of the total correction at T = 4 MeV, its relative 
importance decreases linearly to a 30 % contribution for the T = 20 MeV case. 

Finally, the entropy of Eq. ( 1TB1) is given in terms of the fully correlated momentum 
distributions. This momentum distribution includes both thermal effects (which are capital 
for any entropy computation) plus correlation effects. In fact, in the final entropy S NK both 
effects are taken as the same and correlations somehow mimic extra thermal effects. This 
is why this is the only approximation which tends to give a non-zero entropy at T = 0. 
In the fully degenerate limit, the momentum distributions given by Eq. ( fT9l) are not Fermi 
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step functions and they are corrected by correlation effects. These correlations are the 
responsible for a certain amount of entropy, when this is computed with Eq. f!78|) at T = 0. 
Thus, at finite temperatures thermal effects are overestimated in S NK due to the presence 
of correlations, and S NK produces a far too large entropy (almost a factor of three too large 
at T = 5 MeV). 

After having computed the entropy, we would like to address the subject of TD con- 



sistency. From first principles, the ladder approximation is known to be ^-derivable [51]. 
Indeed, we have computed an expression for the entropy within a formalism that preserves 
$-derivability. Thus, whenever S' is negligible in our approach, we expect our results to pre- 
serve TD consistency. Figure [TU] shows the accuracy that we reach with our SCGF results. 
For the sake of comparison, we also show the BHF results. The upper solid lines with full 
circles correspond to the free energies per particle computed within the SCGF approach: 

F SCGF = E GMK _ TS DQ ? (g5) 

with the energy computed with the GMK sum rule Eq. (1211) and the entropy with the DQ 
expression Eq. (H5|) . The BHF free energy is shown in a full line with solid diamonds and is 
simply given by: 

p BHF = E BHF _ TS BHF ^ (g6) 

with the energy computed from the generalization of the T = BHF approach Eq. fl30|) and 
the entropy from Eq. (1321) . The dotted lines with empty circles and diamonds correspond to 
the microscopic chemical potentials \x obtained from inverting Eq. fTTTl) for the SCGF and 
from inverting: 

P=vj ^f^BHF(k),t BHF ) (87) 

in the case of BHF, respectively. Both chemical potentials are compared with the macro- 
scopic chemical potentials, fi, obtained from the derivatives of the free energies with respect 
to density Eq. (1721) . shown with a dashed line for the SCGF and by a dashed-dotted line for 
SCGF. The derivative has been performed numerically after adjusting F to a third-order 
polynomial. Although the low density region is not well reproduced in this rough approxi- 
mation, the results of the intermediate density region can be fully trusted and in addition 
they are smooth with density. 
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The fulfilment of TD consistency for the SCGF approach is nicely illustrated in Fig. [TU1 
The agreement between p and p is very good above 0.05 fm~ 3 , with discrepancies of less 
than 1 MeV up to p = 0.5 fm~ 3 . As a consequence, the Hugenholtz-van Hove theorem is also 
very well fulfilled, and the minimum of F/A and p do nicely coincide at about p ~ 0.27 fm -3 . 
The situation for the BHF approach, on the other hand, is much worse, as it is very well 
known j^]. The chemical potentials p and p differ by about 10 MeV at p = 0.16 fm~ 3 and 
by almost 30 MeV at the highest density here considered. In addition, the Hugenholtz-van 
Hove theorem is badly violated, and the value of F/A at saturation differs from p by about 
20 MeV. Finally, we note that the propagation of holes seems to have a global repulsive 
effect on the free-energy. Such a repulsive effect has already been detected for the energy 
per particle [33J]. It is interesting to note, however, that the differences in free energies 
between the BHF and the SCGF are smaller than the differences in energies, due to the 
different entropy contributions. 

The results presented here should be taken as a first step towards a full treatment of the 
thermal properties of infinite nuclear matter within a formalism that includes short-range 
correlations in a TD consistent way. From our point of view, this formalism can find several 
applications within the many-body and nuclear physics community. A first outcome could 
be, for instance, the study of the liquid-gas phase transition of symmetric nuclear matter 
from a realistic NN potential. A critical study of the usually assumed low-temperature de- 
pendences for the relevant TD properties of matter (of the type T 2 ) could also be assessed. 
Other applications, such as the study of the finite temperature equation of state, will proba- 
bly demand for an inclusion of three-body forces in the formalism in order to reproduce the 
empirical saturation properties of nuclear matter. As for the possible consequences of our 



results for the study of intermediate energy heavy ion collisions [64 1. we expect that they 
could affect the analysis based on transport models. The use of spectral functions in the 
description of correlated nucleons goes beyond the quasi-particle picture customarily used in 
transport codes. Although some results point out the fact that off-shell effects in the prop- 
agation of nucleons are small 65j, a treatment of a kinetic equation includin g fu ll spectral 
functions obtained from realistic NN potentials (following, for instance, 66j or 67J]) is, to our 
knowledge, still lacking. Moreover, even within the usual quasi-particle description, some of 
the in-medium modifications of nucleons (such as effective masses or NN cross-sections) are 
taken as simple parameterizations [681 ] . Our model permits the calculations of these quan- 
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tities from realistic NN potentials in a fully microscopical and TD consistent basis which, 
properly parameterized, could be used in this kind of studies. 

V. SUMMARY AND CONCLUSIONS 

The ladder or T-matrix approximation of the SCGF method is a $ derivable approach 
and thus theoretically it should fulfil TD consistency. We have checked the numerical TD 
consistency of this approach for symmetric nuclear matter at finite temperature by com- 
puting an approximation to the entropy, the dynamical quasi-particle entropy, S D ®, which 
has been discussed by Carneiro and Pethick Using this S D ® approximation we obtain a 
good agreement between the chemical potentials determined within the SCGF calculation 
and the corresponding values derived from the thermodynamical relations. Therefore, the 
ladder approximation of the SCGF approach supplemented by the evaluation of the entropy 
gDQ provides a method of calculating thermodynamic properties of nuclear matter, which 
accounts for correlations beyond the mean field picture in a consistent way. 

The entropy, S D ® , can be evaluated in terms of a weighting function B, which is connected 
to the usual single-particle spectral function A of the SCGF. This means that B can be 
calculated in term of the nucleon self-energy and an explicit evaluation of the generating 
functional $ can be avoided. The B spectral function and its momentum, energy, density 
and temperature dependences have been studied. In general, one can say that this function 
exhibits a more pronounced quasi-particle structure than the corresponding spectral function 
A. As a consequence, the correlation effects related to the broadening of the quasi-particle 
peak are not very important and a quasi-particle approximation to the evaluation of the 
entropy is a very good approximation if the quasi-particle energies are derived from the 
SCGF self-energy. Even if the quasi-particle energies are approximated by the single-particle 
energies derived from the temperature dependent BHF approximation, the values for the 
calculated entropy deviate only by as much as 10 to 20 percent. In contrast to the SCGF, 
however, the BHF approximation fails to fulfil thermodynamic consistency. The microscopic 
and the thermodynamic chemical potential deviate substantially in the case of BHF and the 
Hugenholtz-van Hove theorem is violated by more than 20 MeV. 
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T [MeV] 


N A (T) [MeV -1 fm -3 ] 


N B (T) [MeV _1 fm- 3 ] 


N QP (T) [MeV^fm- 3 ] 


in 

m 


4 


0.00435 


0.00608 


0.00608 


0.935 


6 


0.00430 


0.00585 


0.00586 


0.901 


8 


0.00424 


0.00566 


0.00570 


0.875 


10 


0.00416 


0.00548 


0.00557 


0.855 



Table I: Densities of states related to the A (first column) and B (second column) spectral functions 
at p = 0.16 fm -3 for different temperatures. The quasi-particle approximation to the density of 
states, Eq. (183 j) . is displayed in the fourth column, together with the effective mass at the Fermi 
surface in the fifth column. 
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Figure 1: Diagrammatical representation of the T-matrix within the ladder approximation. 




I i i f i 
i i i i i i i 



i i i 



i i i 



Figure 2: n-th order contribution to the ladder self-energy. 



35 




Figure 3: functional for the T-matrix approximation. 
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p=0.16fm 3 , T=10MeV 
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Figure 4: (Color online) B (solid lines) and A (dashed lines) spectral functions at p = 0.16 fm 
and T = 10 MeV for three different momenta k = 0, /cf and 2kp. 
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Figure 5: (Color online) B spectral function for five different densities (from p = 0.1 to p = 0.5 
fm~ 3 in equidistant steps) at a fixed temperature of T = 10 MeV and three different momenta 
k = 0,kF and 2Jzf- 
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T=10 MeV 




k [fm" 1 ] 

Figure 6: (Color online) Momentum dependence of the £ function (see Eq. (I75p ) for five different 
densities (from p = 0.1 to p = 0.5 fm~ 3 in equidistant steps) at a fixed temperature of T = 10 
MeV. The dots correspond to the quasi-particle approximation (Q p for the same densities and 
temperature. The arrows signal the position of the Fermi momentum at each density. 
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p=0.16fm 




oo-|i [MeV] 



Figure 7: (Color online) B spectral function for five different temperatures (from T = 4 to T 
MeV in equidistant steps) at a fixed density of p = 0.16 fm -3 and three different momenta k = 
and 2k f- 
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Figure 8: (Color online) Different approximations to the entropy as a function of density for a 
T = 10 MeV temperature. The full lines correspond to S D ®; the dotted lines to S® p ; the dashed 
lines to S BHF and the dot-dashed lines to Si 4 . 
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T [MeV] 



Figure 9: (Color online) Different approximations to the entropy as a function of temperature at 
a density of p = 0.16 fm~ 3 . The full lines correspond to S D ®; the dotted lines to S® p ; the dashed 
lines to S BHF ; the dot-dashed lines to and the double dot-dashed lines to S NK . 
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Figure 10: (Color online) Free energies per particle (full lines) and fx chemical potentials (dotted 
lines) of the SCGF (circles) and BHF (diamonds) approaches as a function of density at a tem- 
perature of T = 10 MeV. The \i chemical potential obtained through a numerical derivative are 
displayed in a dashed line for the SCGF results and a dot-dashed line for the BHF results. 
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